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1. INTRODUCTIOW 


In 1973, the Health Applications Group (HAG) at the National 
Aeronautics and Space Administration (NASA) Johnson Space Center 
(JSC) began a program to evaluate the use o£ remotely sensed 
data as an additional tool in existing and projected insect 
eradication efforts. The need to relate remotely sensed data to 
screworm infestations resulted in a two-part developmental 
study. First, the ability to use remotely sensed data to 
estimate weather conditions would be evaluated. Second, the 
effect of weather on screwworm populations would be modeled. 

The first study, concerning the ability to use remotely sensed 
data, resulted in several reports indicating that the Improved 
TIROS Operational Satelllite (ITOS) appeared to be the most 
desirable remote sensing satellite and that the Johnson Space Cen 
ter had the necessary capability' to process and register the data 
The second study, concerning the screwworm and the role of the 
weather in the activity and development of the screwworm fly, 
also resulted in several reports. Information and illustrations 
from these reports will be included in the following text. Con- 
currently with the developmental studies, a production processor 
was designed and built. 

1.1. BACKGROUND 

The initiative for the joint NASA-U.S. Department of Agriculture 
(USDA) project on the eradication of screwworms was derived from 
Dr. Charles Barnes, PH. D, D.'V.M., of NASA Johnson Space 
Center's Health Applications Group, Life Sciences Directorate. 

He was intensively familiar with the screwworm problem, having 
served earlier as an inspector with the USDA in Mexico. In 
addition, he was familiar with the untapped potential of using 
satellite data for health applications work, such as in support 
of screwworm eradication. 
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The Lockheed Electronics Company, Inc. (LEG)., a support contrac- 
tor, hired a group o£ multidisciplined scientists to provide 
space technology techniques to support an eradication program. 
Processing techniques were developed and test sites were estab- 
lished in Veracruz, Mexico and in Texas for purposes of 
collecting ground truth data. A production processing system, 
Screwworm Eradication Data System (SEDS) , was designed and built 
by the NASA JSC Ground Data System Division (GDSD) and their 
support contractor Aeronutronic Ford. The basic algorithms were 
developed by the HAG/LEC and incorporated with some modifica- 
tions into the SEDS design by GDSD/Aeronutronic Ford. Due to 
concern for operational problems, there was little active 
participation in the project by the USDA until after the initial 
algorithms were developed and the production system was com- 
pleted and ready for evaluation. 

1.2 PURPOSE 

This report documents the development of the basic algorithms 
and explains how the various data are prepared for input into 
SEDS, as well as how operational coefficients are calculated. 
This report provides a knowledgeable programmer the necessary 
and available information with which to build a stand-alone 
system. It is intended to be used in conjunction with the 
basic SEDS documentation for purposes of technology transfer. 

1.3 DEVELOPMENT OF THE SYSTEM 

The development of the SEDS system was limited by the lack of 
reliable biological data for modeling. It was also limited by 
lack of funds, time, and availability of major pieces of compu- 
ter equipment at JSC. Probably the most limiting factor in the 
SEDS processing system was the limited amount of computer 
memory, a factor which hampered the number of functions SEDS 
could perform. 
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The development of SEDS was and still is a cooperative effort; 
the operation of SEDS is also cooperative in that Aeronutronic 
Ford operates the SEDS on a production basis and LEG furnishes 
certain periodic operating coefficients and background correc- 
tion maps which are produced for SEDS on other computers. It 
must be stressed that SEDS could be a stand-alone system provi- 
ded that it had sufficient peripherals, memory core., and suffi- 
cient time for programming these other functions into SEDS. 

1.3.1 DATA AVAILABLE FOR ANALYSIS 

To quantify the effects of a parameter affecting screwworm growth, 
a statistical relationship between that parameter and the screw- 
worm population must be shown. Among the important parameters 
that affect scre'wworm population dynamics are the availability 
of hosts, incidence of wounds, predators, efforts of the eradi- 
cation program, daily temperatures, and soil moisture budget. 

Of these parameters, data on only the latter two were available 
and statistically reliable enough to warrant detailed analysis 
and comparison with screwworm population changes. Although 
screwworm populations could not be determined, the population 
changes were inferred from the number of screwworm cases. 

The following data was used in the development of the screwworm 
growth potential functions. 


Year 

State 

Meteorological Divisions 

1968 

Texas 

5,6,7,9,10 

1973 

Arizona 

1,3, 4, 6, ,7 

1973 

Texas 

5,6,7,9,10 

1974 

(to 31 July) 

Texas 

5,9 


Data was not analyzed outside of Texas except for Arizona (1973) . 
The number of screwworm cases reported in 1967, 1969, 1970, and 
1971 was low and was not used. The 1972 data represented a 


1-3 



large anomalous outbreak that remains largely unexplained, and 
was not used to develop the growth potential functions. A 
later test of Division 9 (1972) revealed that at least that di- 
vision could have been used without affecting the final results. 
The data was censored to ten cases qr more. Below ten cases, 
the data was not considered reliable, and therefore, was not 
used . 

The meteorological data consisted of weekly mean air tempera- 
tures, total precipitation and crop moisture indexes for each of 
the divisions. The values used represent the average of data 
from all of the reporting stations within each division. 

The screwworm cases were converted to growth per generation, 
that fs , growth is defined as the ratio of the number of cases 
of a given week to the number of cases reported during the \^eek 
in which the parent generation occurred. 

A weekly plot was made of the meteorological data versus growth 
per generation as shown on figure 1-1. A curve was then hand 
drawn in a manner which, in the analyst’s opinion, most closely 
fit the relationship of the two parameters being plotted. Each 
curve defined the growth potential for that meteorological para- 
meter. Growth .potential of a parameter is defined as the maximum 
likely growth, if that parameter is the only factor restricting 
screwworm growth. 

1.3.2 DATA USED TO DEVELOP SCREWWORM MODEL 

Extensive literature reviews of the effect of environmental con- 
ditions on the growth and development of the screvrworm were used 
as the point of departure in producing a screwworm model. A 
brief review and summary of the surveys is contained in the 
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technical memorandum, "Selected Parameters Influencing the 
Growth Survival of Cocht'iomyi.a honvln'ivorax (Coquerel)." 

The generalized life cycle of the screwworm fly is presented in 
figure 1-2. Each stage of the life cycle has key factors which 
influences its development rate and which would have been desir- 
able to incorporate into the model. However, the only statisti- 
cally reliable data was on weather. A more complete description 
of the effects of environmental factors at each stage and some 
consideration for modeling is contained in the above-referenced 
technical memorandum. 

The relationship of crop moisture index (CMlj and dally mean air 
temperature (DMAT) is clearly illustrated in figure 1-3. It 
shows that both temperature and CMI are the primary limiting 
factors affecting screwworm growth at different times of the 
season. One of the modeling considerations is the fact that the 
daily air temperature affects both the adult fly and the pupae 
in the ground. The newly hatched or adult fly are immediately 
affected by the existing air temperature, whereas the pupae are 
not, or only slightly, affected by the then existing air tem- 
perature. Pupae are affected more so by the daily mean air tem- 
perature over an extended period due to their slow development. 
Therefore, specifications for SEDS included multiple output 
products: long-term growth potential, short-term growth 

potential, CMI and degree day summation. These output products 
were to be used by trained USDA epidemiologists to assist them 
in planning fly drop missions. The epidemiologists were to have 
had the detailed knowledge of their area and other factors 
affecting screwworm growth, so that they could intelligently use 
the additional tools depicting environmental conditions espe- 
cially in remote areas. 
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Figure 1-3.— Data plot showing the number o£ 
screwworm cases and the weather in Texas 
Climatological Division 9, 1968. 






At a later stage in the development of SEDS, a decision was made 
that the four output products should be combined. A regression 
analysis of the four products was performed, and a preliminary 
set of coefficients was calculated for combining the four -pro- 
ducts. It is expected that these coefficients will change for 
different seasons of the year. It may also be possible after 
some experience to eliminate one or more output products without 
a loss of overall results. 

1 . 4 EVALUATION OF THE SYSTEM 

SEDS first produced data in March 1975, and there was an immedi- 
ate attempt to evaluate its effectiveness. During this initial 
evaluation, it became readily apparent that there needed to be a 
temperature fill-in scheme incorporated into SEDS to provide 
data when clouds obscured the ground. It was decided, there- 
fore, to incorporate a delta temperature field and other fixes 
into SEDS and reprocess the spring 1975 data. An evaluation 
team of NASA and USDA representatives began their evaluation on 
the data as it was produced weekly. 

Completion of spring data reprocessing occurred on January 9, 
1976. The ongoing evaluation was accelerated, so that most of . 
it could be reported at the following meeting of the full eval- 
uation team. 

The results of the spring evaluation are contained in separate 
documents or reports and are beyond the scope of this document. 
There were two factors associated with the evaluation which are 
worthy of listing since they could apply to any future evalua- 
tions. The first major limiting factor in the evaluation was 
the relatively small number of cases in the early spring from 
which to calculate actual growth potential for use in statisti- 
cally comparing with SEDS growth potentials. For this reason, 
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much of the evaluation had to be subjective. It would have been 
preferable to have had sufficient data points from which to do 
a statistical evaluation. 

The evaluation of the reprocessed spring data was also adversely 
affected by the speed at which it was produced. Since the 
registration coefficients had already been previously recalcu- 
lated, that portion of the data did not need to be repeated. 
Consequently, four days of spring data could be reprocessed in 
each production day, leaving insufficient time to adequately 
analyze the data and compute the necessary operating coeffici- 
ents for SEDS to perform to its capacity during the spring 
evaluation period. 
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2, CALCULATION OF SEDS SYSTEM PARAMETERS 

The daily operation o£ SEDS requires the availability of several 
sets of parameters. These parameters are input to the system 
from a list of currently valid values. The SEDS system para- 
meters form shown in figure 2-1 is also used to initiate any 
changes in system coefficients. 

The following sections provide details on the calculation and 
utilization of each set of coefficients. Whenever possible, 
information has been provided on the theory and development of 
the current configuration of SEDS. However, emphasis is on 
developing a complete understanding of the analysis procedure 
required to maintain the highest possible system performance. 
Examples are provided both for calculation of the coefficient 
sets and for determining when recalculation is needed. 

2.1 DAILY MEAN AIR TEMPERATURE (DMAT) COEFFICIENTS 

SEDS uses a set of coefficients derived from multiple linear 
regression studies to convert radiometric data obtained by the 
National Oceanic And Atmospheric Administration fNOAA) satellite 
into estimates of mean air temperature. Computed daily, the 
mean air temperature is used by the various SEDS screwworm 
population dynamics algprithms to estimate the effect of tem- 
perature on different phases of the screwworm’ s life cycle. The 
performance of the SEDS screwworm growth potential predictions 
depends largely upon the accuracy of the estimated daily mean 
air temperature. 

2.1.1 BACKGROUND MATERIAL 

The DMAT estimate is calculated by SEDS each day for each pixel 
in the data base image. The exact method of calculation depends 
upon the amount of radiometric data available. Data 


2-1 



SEDS SYSTEM PARAMETERS 


Disk Number 


Date 


DMAT Coefficients 



10 - 


^ 20 - 


" 30 - 


" 11 - 


" 21 - 


" 31 - 


" 12 - 


" 22 - 


"32- 


" 40 - 


"50- 


" 41 - 


" 51 - 


" 42 - 


Day Thresh 


" 52 - 


Night Thresh 


Screwworm Growth Coefficients 



Ground Truth Constants 

KDAY MINDLT 

KNIGHT DDT 

Isothermal Bias and Offset 

A B 

*Parameters in parenthesis indicate those which are not in converted 
(internal) units. 

Programming Support Section 

Earth Resources Operations Section 


Figure 2-1. - SEDS System Parameter Control Form. 
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availability is largely a function of satellite coverage and of 
cloud cover. 

The radiometric data must survive a number of filters which 
strive to eliminate all cloud contaminated data. Although 
these filters are described elsewhere, it is worth emphasizing 
that the quality both of the coefficients determination and of 
the DMAT estimates is dependent upon the use of cloud-free 
data. The quality of the registration of the radiometric data 
to the SEDS reference image is also extremely important. 

Without accurate registration, no meaningful comparison can be 
made between surface mean air temperature and the satellite 
observed radiometric temperatures. 

2. 1.1.1 Theory and Development 

Two assumptions are involved in attempting to relate radiometric 
temperatures to daily mean air temperatures . The first is that 
the instantaneous radiometric temperature, which is a surface 
temperature, is related to the air temperature measured at some 
fixed height above the ground. The second assumption is that 
one or two estimates of air temperature each day will be suffi- 
cient data to estimate the daily mean air temperature. Although 
in reality neither assumption is completely valid; as approxima- 
tions of the complex, real situation, they appear to provide 
usable data. 

On an experimental level, it can be seen from figure 2-2 that 
in data tahen from a test area in Veracruz, Mexico, there was a 
good (.811^ correlation between radiometric and air temperatures 
taken at the time of satellite passage. Simulations using sur- 
face observed air temperatures at the time of satellite passage, 
indicated that the potential for predicting daily mean air tem- 
peratures was quite good. Based on these data from ten Texas 
stations from the early summer of 1973, it was found that by 
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SURFACE OBSERVED AIR TEMPERATURES (®C) 



SATELLITE OBSERVED RADIOMETRIC TEMPERATURE (°C) 

jFigure 2-2. - Comparison of satellite data with surface air 
temperature in the Fortin de las Flores test area. 


2-4 



knowing the air temperature perfectly the daily mean could be 
predicted with error ranging from 1.07 to 3,43°C, depending upon 
the availability of data. 

It can be shown that the effect of the air -ground temperature 
difference dominates the error produced by estimating actual 
surface net radiation flux from data taken in a weather shelter 
above the ground. In this case, the actual surface temperature, 
which is a function of the surface net radiation flux, is 
measured and an estimate is made of the net radiation flux (ur 
temperature) of the air. If the air and surface are in thermal 
equilibrium, the temperature difference remains much the same 
over relatively short periods of time and under the influence of 
similar weather conditions. 

Thus, it should be expected that radiometric data taken at night 
under relatively stable conditions would provide more reliable 
air temperature estimates than data taken in the early morning 
when the surface is undergoing rapid solar heating. Addition- 
ally, since the air-ground temperature difference has a strong 
seasonal dependence, it should be expected that any transforma- 
tion between radiometric and air temperature would change with 
time. 

A series of regression studies performed on July 1975 SEDS data 
revealed that the addition of terrain elevation as a variable 
could improve the errors in estimating mean air temperature on 
the order of 1°C. Tabl-e 2-1 shows the relative improvement 
gained by including the additional variable. The data agrees 
with the expected results quite well. The DMAT estimates based 
solely on day radiometric data are of a lower quality than those 
based on night data. The elevation correction has little effect 
on the errors associated with the high quality night data. 
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TABLE 2 - 1 , - IMPROVEMENT OF DMAT ESTIMATES UPON ADDITION OF 
ELEVATION AS A REGRESSION VARIABLE TO JULY 1975 SEDS DATA 



Radiometric 

Altitude 


Data Only 

Corrected 

Day Radiometric Only 

6.5 

3.4- 

Night Radiometric Only 

2.6 

2.0 

Both Day and Night Radiometric Data 

3.6 

2.1 ■ 

Overall (All Radiometric Data) 

3.48 

2.39 


2. 1.1. 2 SEDS Utilization 

Depending upon the availability of radiometric data, there are 
three possible cases as shown below. Usable data is radiometric 
data which has survived all cloud detection algorithms. 


» Case 1: Both night and day radiometric data usable 

DMATK = A + A-TSNK + A„TSDK + A„ALT (2-1) 

0 1 ct o 


9 Case 2: Only day radiometric data usable 

DMATK = B + B.TSDK + B„ALT 
0 1 2 

• Case 3: Only night radiometric data usable 

DMATK = + C^TSNK + C^ALT 


( 2 - 2 ) 


(2-3) 


Where the symbols are defined as follows ; 

DMATK — Daily mean air temperature in degrees Kelvin 
TSDK — Day radiometric temperature in degrees Kelvin 
TSNK — Night radiometric temperature in degrees Kelvin 
ALT — In meters (terrain elevation) 

A., B,, C. — Regression coefficients. 

'h' "tf 
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Computationally SEDS uses the following relationships: 
e Case 1: DMAT = AOC + AICCTSN) + A2GCTSD) + A3C(PALT) (2-4) 

• Case 2: DMAT = BOC + BIC(TSD) + B2C(PALT) (2-5) 

• Case 3: DMAT = COC + CIC(TSN) + C2C(PALT) (2-6) 

Where the following symbols are used: 

• DMAT — Eight-bit coded value of DMATK. DMATK and DMAT are 

related by the following 

DMATK = DMAT + 260 
~T~ 

9 TSD — Eight-bit coded value of TSDK. TSDK and TSD are 
related by the. following 

TSDK = TSD + 202 
2 

• TSN — Eight-bit coded value of TSNK. TSNK and TSN are 

related by the following 

TSNK = TSN + 202 
2 

0 PALT — Eight-bit coded value of ALT. ALT and PALT are 
related by the following 

ALT = 25 (PALT) - 100. 

The coded coefficients are floating point numbers which satisfy 


the following: 

AOC = 4A + 808A- + 808A„ - 400A„ - 1040 (2-7) 

O -L a O 

AlC = 2A^ (2-8) 

A2C = ZAg (2-9) 

A3C = lOOAj (2-10) 
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BOC = 4B^ + 808B^ - 400Bg ~ 1040 (2-113 
BIG = 2B^ (2-12) 
B2C = lOOB^ (2-13) 
COG = 4G + 808G- - 400G„ - 1040 (2-14) 
CIC = 2C^ (2-15) 
C2G = lOOC^ (2-16) 


2.1.2 METHOD OF CALCULATION 

The general approach to calculation is to develop a data set for 
each of the three cases. The data sets should contain cloud-free 
radiometric 'data which cover the entire range of expected values. 
Multiple linear regression analysis is used to calculate the 
coefficients for equations (2-1), (2-2), and (2-3). The coeffi- 
cients are transformed into the form- expected by SEDS and entered 
into the systems parameter list for implementation. A given set 
of coefficients remains valid for 4-6 weeks. In general, a 
continuing error analysis indicates when recalculation is 
necessary . 

2. 1.2.1 Data Used 

The radiometric data source is the SEDS meteorological station 
statistics printout. A typical printout is shown in table 2- II. 
The values for TSDK, TSNK, and VIS for those stations in the SEDS 
prime verification network are extracted. The stations in the 
prime verification network are shown in table 2- III. In general, 
at least two calendar weeks of SEDS data are required to be sure 
of obtaining a sufficient amount of data. 

The primary source of surface meteorological data is the daily 
report provided by the National Weather Service. An example of 
this report is shown in table 2- IV. The daily mean air 
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MGINAi; PAGE IS 
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TABLE 2-II. - PARTIAL LIST OF METEOROLOeiCAL STATION STATISTICS 

PRINTOUT FOR 31 MARCH 1975 


ID 

Y 

X 

TSD 

TSN 

DMAT 

R 

CMI 

DAY TMET 

DELT TGT 

ALT 

FL/S 

CODE GRID 

GRID 

K 

K 

K 

CODE 

CODE 

VIS K 

K 

K 

M 



« 

1U 

251,0 

0.0 

2*8‘.25 

M 

0,18 

118 

8*.7,00 

1,8$ 

<*•.83 

750,0 

n 


12 

38a 

0*7;5 

3, A 

2*7' 75 

CO 

A. 18 

44 

2*6,60 

1.83 

8*7,7$ 

130,0 

n 

at. 

- 11 

158„ 

-Jbt.S 

— 6.0 

290.00 . 


--0.U - 

. 49 

8*7,00 . 

-- 1,00 

890.00 

. 430,0 

a 

«T 

19 

3*5 

877.0 

0,0 

1 * 8.75 

k 

0.12 

13 

2*6,50 

8,8$ 

8*8.7$ 

130,0 

n 


10 

2*4 

2*7. » 

0.0 

2 * 6.25 

CO 

6.18 

• 9 

8*6.60 

>4.8$ 

2«6.8$ 

T 56 .O 

0 


u 

?5» 

i9*:o 

0.0 

2*6.75 

CD 

8.18 

«1 

2*5,30 

1 . 8 $ 

8 * 6 . 7 $ 

730.0 

0 

12 

It 

84) 

2Ta;5 

0.0 

2*2.75 

K 

0.18 

U 1 

8*5.2$ 

•0 ,30 

282.7$ 

150.0 

0 


19 

482 

2*2,5 

n.n 

2*1.75 

K 

l .*0 

ISO 

2*5.8$ 

0,50 

2 * 5 . 7 $ 

150,0 

0 

1 « 

-.42 

500 

-864,$ - 

- 0,6 - 

8A7;$0 

-k- - 

--o.ia 

-$5 

8*4,60- 

.— - 1,00 

887,50 

760,0- 

-.a 

1 * 

42 

48A 

2«2.5 

0,0 

2 * 1.00 

M 

!.*0 

1 5* 

2 * 1 . 2 $ 

- 0 , 2 $ 

8*5.00 

130 .0 

0 


t* 

207 

29*. 5 

0.0 

2*6.60 

CO 

•e, 18 

44 

8*6,50 

1.00 

886,30 

1285,0 

a 


49 

447 

|7*.5 

0.6 

2*2,75 

■ H 

0,12 

117 

2*5.8$ 

-0,50 

8*2.7$ 

130.8 

0 


15 

559 

291,5 

0,0 

2*7,59 

* 

0,18 

19 

8*6.30 

8.50 

289,8$ 

450,0 

0 

^TJ 

IS 

105 

|90.0 

0.0 

2*4. 0(1 r 

CO 

«0.3T 

48 

8*8.50 

1,30 

2*4,00 

1885,0 

0 


$4. . 

. .at t. - ^ - 

. a a 

2*6',25‘ 

'.L 

A t?— 

- jkT aAT- Afl_ 

- 0,75 - 

- jSa, 

l$0_ft 

_o 

7a 

b2 

142 

—9# 

2«3.5 

6,0 

2*6,60 

CO 

— T 

-0,37 

6$ 

2*5,50 

1 .00 

8*6.60 

128$. " 

a 

... S« 

. bS 

lat . - 

29$. 0 

0,0 

286.50 

CO 

.0.37 

52 

2*5.30 

1.00 

<86.50 

128$. 0 

0 

«0 

k5 

41* 

267.$ 

0.0 

2*7,00 

N 

A, 12 

112 

8*7,00 

0 . 0 c 

8*7.00 

150,0 

a 

- ,86 

b7 

2S7 . 

297,0 

0,0 

2«b,75 

CO 

-0.7S 

«t 

8*5.50 

1.25 

2*6.75 

1283,0 

0 

-TO , 

Ta 

21 

0,0 

0.* 

2*7.25 

3 

0.12 

0 

2*».50 

0.75 

2*7.2$ 

750.0 

0 

A 

_ 7b 

.178 

-iSI.O-- 

0,0 . 

-2*7,00 - 

CO 

_»0,3T_ 

—UB 

-2B5,$0 

1,50 

247,00 

1225,0 . 

..0 


79 

429 

274.$ 

6.0 

2*8,75 

H 

l.*!0 

1«3 

8*0.2$ 

0.30 

8*4,75 

ISO.O 

0 

aa — 

85 

29S 

298,0 

0,0 

288.75 

CO 

-1,85 

43 

2*9.2$ 

-0 .50 

288,75 

450,0 

0 

83 

8a 

S6* 

2Tl'.$ 

0.0 

2*7.75 

L 

0,12 

49 

8*7.00 

0.75 

8*7,75 

150,0 

a 

.. .ba . . 

87 

389 

274', 8 

0,0 

28r'.50 

H 

.0,12 

183 

8*4.85 

-0.75 

8*7,50 

150.0 

0 

Tl 

B« 

39 

0.0 

0.0 

2*S;00 

0 

0.12 

0 

2*1.0* 

2.00 

8*5.0" 

1225.0 

0 

_.^J2 

_89_ 

_ 60 — _ 

.291,0-- 

--0,0 . 

-2&3.2S 

CO- 

-«0,12 - 

-45. 

-8*2.00 . 

- 1.25 

885,2$ 

.5225,0 

. a . 

M 

90 

259 

292.5 

0,6 

2*8.5" 

CO 

-1.85 

ag 

2*0.6" 

8,00 

2*6,50 

760,0 

0 

_ — 75 _ 

. 95 

215 

29710 

0,0 

2>‘7.75 

c 

-0.37 

35 

2*0.50 

5.00 

8!6.2$ 

1285.0 

0 

So ' 

«a 

237 • 

294.$ 

i.9 

2*5.25 

CD 

-0.75 

ao 

2*0.50 

4,7$ 

8*3.25 

1285.0 

0 

.. '^’T 

«a 

«IS . 

254.0 

0.0 

248’. So 

H 

.0.12 

143 

2*4.2$ 

0.25 

2*4.50 

iso.fl 

0 

55 

102 

24 

0,0 

6.0 

2*7.75 

0 

0,12 

0 

2*6.50 

1.2$ 

287.7$ 

bOO.O 

0 

'^76 

.105_ 

304 

_22fllo 

_ 2fl9'^25 

CO 

- -1 .25 

—90 

_2B9,25_ 

O.OQ - 

269.25— 

- a5fl..o„ 

D- 

e« 

112 

397 

279,0 

6,0 

249;$0 

H 

-0.75 

140 

2*8,25 

1.25 

249.5" 

150.0 

0 

15 _ - 

520 

210 

300.5 . 

0.0 

288.00 

CO- 

-0,37 . 

. 4b . 

. 2*0.50 

.. -7.50 

2*8.00 

1885.0 

0 

--77 

1?2 

356 

282.0 

0,0 

248.50 

CO 

-0.37 

51 

2*6,2$ 

0.25 

8*8.50 

ISO.O 

0 

-. ?7 . 

.125. 

.291. .. 

,29 7',0 . 

. 0,0 

288’.75 

■C _ . 

.-0.75 

39 

2*0. SO 

S.SO 

8*5.5" 

7SO.O 

0 

5b 

125 

121 

299.5 

e.A 

291,25 

c 

-0,75 

33 

2*9.25 

t,50 

890,75 

150,0 

0 

-«b 

-15 1-. 

SO* 

.£14*5—. 

fl,6 

2?Ql50- 

' u 

-0 , 75 

_99_ 

.£*9..2S_ 

I-J5.- 

.290^51) _ 

— 1.511. 0_ 

__ n.- 

B9 

152 

”357 

302.5 

0.0 

292,75 

c 

-0,75 

58 

288.25 

0.50 

288.25 

150.0 

0 

.. -.52 

555 

356 

300.0 ■ 

..0.0 

291.50 _ 

.CL_.. 

-0,75 

37 

2*8.25 

0.25 

2*6,00 

150.0 

0 

80 

153 

380 

100.0 

0.0 

291,56 

c 

-0,75 

33 

2*9.25 

1 .50 

290 .SO 

150.0 

0 

. Z ^ 

.158 

171 

297.5 . 

0,0 

288.50 

CD 

-5.37 

. 52 

2*9.25 

-0.75 

2*8.50 

8500.0 

0 


156 

408 

246,5 

0.3 

290^53 

M 

-l.SS 

ll9 

2*9.2$ 

1.25 

290,50 

ISO.O 

0 

^ 1 

.169_ 


^ fl 

0 J1-. 

29 L»0t 

_6_. 

_.0-12 

0. 

290, 2S_ 

0,-15.. 

_2?1.0'i.. 

baO.Q,— 

_B.. 


179 

'4I8 

246.5 

0.0 

2* 8'. 25 

K 

-1 .25 

12* 

8*6.75 

-0.50 

2*8,25 

ISO.O 

0 

. — 78 

.18a 

400 . 

2Sa;o 

0,0 

2*6.75 

M 

. .1.25, - 

105 . 

2*9.75 

- -1 ,00 

288.7.5 . 

. 150,0 - 

-0. 

79 

187 

4^4 

257.0 

n.o 

2*9,03 

H 

-1,25 

13* 

2*8.7$ 

0.25 

249 ,00 

150,0 

0 

87 

186 

423 

.256.0 

0,0 

249.00 . 

... H 

-1,25 

142 

2*8.75 

0.25 

2*9.00 

150,0 

0 

57 

189 

40? 

24 b;. 0 

0.0 

286.25 

M 

-1.26 

126 

2»9.7$ 

-1 .50 

2*6.8$ 

150.0 

0 

7 


19 

o', n _ 

3,Q.._ 

.28 9-25- 


1_S4— 

^ Q 

-269,2$ 

Q,00_ 

-259.^S_ 

A , 0 

A 

b 

287 

34* 

302,5 

0,0 

292.00 



c 

-0.75 

40 

2*8,25 

0.25 

2*4.00 

bOO.O 

0 

.38 .. 

2l3. 

72 

297,0 

0,0 

.258,75 

.CO. . 

-.0,37 _ 

. «1. 

2*9.25 

• 0.50 

2*8.75 

100. Q 

0 . 

59 

229 

9 

0.0 

0,0 

239,75 

0 

1,50 

0 

2*9.25 

0,50 

2*9,75 

0.0 

0 

-5 

2Sa 

262 

301.5 

0.0 

291'. 7$ 

-CO 

-0,75 

47 

290.25 

1 ,50 

291.7$ 

1500,0 

0 

— 12 

257- 

42l 

246'. 0 

0.0 

2*9.25 

k 

-t.25 

75 

292,2$ 

-3.00 

2*9.25 

100.0 

0 

- 11 

2kO 

194 

290.0 

0.0 

2*5,25 

H 

-1.25 

10b 

290.2$ 

-5.00 

2*5.8$ 

600.0 

0 

> 

282 

157 

296'. 3 

0.0 

293'. 25 

CO 

-*.57 

41 

2*2.25 

1.00 

293.85 

600.0 

0 

-- 9 . 

291 

23* 

300.5 

3.0 

292.75 

CO . 

-0.37 

47 

2*0.25 

.2.50 

298.75 

bOO.O 

0 . 

— 1« 

35 1 

44? 

29$. 0 

0.0 

290.53- 

'CD 

-1.25 

47 

2°2.25 

-1.75 

890,50 

IQO.O 

0 

. bi 

319 

4 1 4 

259.5 . 

0.0. 

292.75 

1 . 

, -1,25 

75. 

294.25 

-1.30 

292, >$ 

IDO.O . 

0 

— 7 

117 

78 

297. S 

0.0 

2*8.53- 

CD 

-0.37 

40 

2*9.2$ 

-0,75 

2*6 .5.0 

100.0 

0 

SB 

J?5 

355 . _ 

ioa.5 .. 

. Jl.3- 

291.54- 

CD— 

— -0..75-... 

- 47.. 

294.25- 

— »2,25. 

-291.50- 

-1.500-Jl _ 


17 

3?5 

St4 

304.0 

0.0 

29l‘ 00- 


.0.75 

37 

29«.2S 

-3.2$ 

2 9 1 0 ^ 

iSoo.o 

a 

- IS 

326 

I9u 

301.5 

. 3,0 

292.50 

c 

.0,37 

33 

292.25 

O.OC 

292.85 

100,0 

e 

- 20 

336 

Ub* 

290.0 

0.3 

290.50. 

H 

.6.75 

iO*i 

292.2$ 

-1.75 

290.30 

100.0 

0 

S6 

342 

319 

29?; 5 

9.0 

291.00 

CD 

.3.75 

5! 

294.25 

• 3.25 

291,0" 

600,0- 

0 

51 

355 

46B 

290,0 

0,0 

290.75 

K 

-0,75 

U? 

292,25 

-1,60 

290,7* 

0,0 

a 

lb 

.381 

34? . . 

.311.5 _ 

5,0 

296. 75_. 

..c . 

. -9,7s 

35 

296.25 

..-2,50 . 

.293.2$ - 

. 600.0 . 

-„0- 

... 15 

589 

298 

285.0 

0.6 

292'. 00 

k 

-0.37 

5? 

292.25 

-0.25 

298. , - 

600,0 

a 

bS 

389 

4U4 

5fta’,$ 

9.0 

2*9, So 

c 

-0,7$ 

U? 

0.00 

• 2.30 

0, 

8500,0 

a 

- 19 

370 

"33 

29*, 0 

3.0 

2*6 00 - 

. C — 

-0.75- 

■ if>: 

0.00- 

-2.00 

0.0" 

2500,0 

0 

85 

37b 

490 

253.5 

0.0 

291,7$ 

K 

.1 .25 

9? 

297.25 

-5.50 

291 .75 

600,0 

0 

— ?b 

i'O 

518 

2*1,5 

0,3 

294.2$ 

CO 

.9,75 

70 

2*7.2$ 

-3,00 

294.8$ 

8,0 

0 
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TABLE III. - SEDS PRIME VERIFICATION NETWORK 


Station Name 

SEDS ID ■ 

Latitude 

'Longitude 

Hermosillo 

01 

29:04 

110:57 

Chihuahua 

02 

28:38 

106:04 

Guaymas 

03 

27:55 

110:53 

La Paz 

07 

29:04 

110:20 

Durango 

09 

24:01 

104:40 

Mazatlan 

10 

23:12 

■106:25 

Soto la Marina 

12 

23:46 

98:12 

Tampico 

14 

22:13 

97:51 

Tepic 

15 

21:30 

104:53 

Rio Verde 

17 

21:56 

100:00 

Guadalajara 

18 

20:36 

103:23 

Pachuca 

19 

20:08 

98:44 

Tuxpan 

20 

20:57 

97:23 

Puebla 

25 

19:03 

98:12 

Veracruz 

26 

19:09 

96:07 

Oaxaca 

30 

17:03 

96:43 

Acapulco 

31 

16:50 

99:55 

Salina Cruz 

32 

16:10 

95:12 

Tucson 

70 

32:07 

110:56 

Dougl as 

72 

31:27 

109:36 

Deming 

73 

32:15 

107:42 

El Paso 

74 

31:48 

106:24 

Marfa 

75 

30:22 

104:01 

Del Rio 

76 

29:22 

100:55- 

Cotulla 

77 

28:27 

99:13 , 

McAl 1 en 

78 

26:11 

98:14 

Brownsville 

79 

25.: 54 

97:26 

Midland-Odessa 

83 

31 :57 

102:11 

San Angelo 

84 

31 :22 

100:30 

Abilene 

85 

32:25 

99:41 

San Antonio 

93 

29:32 

98:28 

Austin 

94' 

30:18 

97:42 

Waco 

95 

31:37 

97:13 

Corpus Christi 

96 

27:46 

97:30 

Victoria 

97 

28:51 

96:55 

Houston 

98 

29:58 

95:12 

Port Arthur 

99 

29:57 

94:01 
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table 2- IV. - DAILY SURFACE DATA REPORT 


^W4vHeXJCO ,;SC|CewWORM~ DATA-P0R-^-:S(;jf^*.2^T^Cvt 1-975 ■ 


^T*f 

• --TT 

RRRR 

ID STN ^ 

:• TT- ' RRRR 

ID STN 

TT 






• 


■M)<ican 

S,T-^T^iCDN:£>AT-A 

FOR ib/27 0753 



oa 122 



•’.M2 499 . 

•^ ;-Ul if ;> 10027 - 

24 680 

18 

01 160 


9'999 

. ' 13. '51 8<; 

'#..:%999- 

25 685' 

■ 17 

02 225 

18 • 

. 0000 

14 -548 ' 

‘0000 

26 692 

28 

03 255 

22 

0000 

15 556 

22' 'OOOO 

27 727 

99 

OA 342 

99 

0000 

16 57-5' 

21 9999 

28 741 

27 

05 382 

.19 

9999 

17 581 

,19 0039 

29 762 ■ 

22 

06 393 

19 

9999 

18 6.13 

0000 

30 775 

20 

07 405 

23 

0000 

1.9 632 

15 0000 

31 805 

27 

08 412 

26 

9999 

20 .640 

19' 0118 

32 833 

27 

09 423 

19 

9999 

21 654 

9.9 ■ 9999 

33 845 

16 

10 -458 

25 

00.00 

22 662 

,99- 9999'. 

34 855 

99 

11. 491 

16 

0000 

23 679 

. 18 • 0000 

. • - 
, . ? 

35 903 

26 

36 AGS 

n . ■ 


48-MLar 

•'.-^99 

60 TCG 

14 

•37 AGU 

■ .99- 


49’ MLM-- 

20 

61 TGZ 

25 

- 3S CEW .. 

21-. 

• 

50 MTT 

- ■. 29 

62 TMN 

17 

39 CJS 



51 NAU • 

25 

63 UPN 

17 

40 GRN 

-.-■99; . 


52 NLD . 

22 

64 VCT 

99 

41 IMA 

, 22 


53, NQG. ’ 

.22 

65 VSA 

27 

42 IZT 

-28.*. 


. ..S4.'FTd. 

. 99 - 

66 ZIH . 

27 

43 JAL 

18 


55 ' FVR- 

^ 27 

80 TIJ 

16 

44 LMM 

99 


56 PZAi, 

.-99- - 

81 MXL 

24 

45 LTD 

99 


57 REX 

^'99-. . 

82 PPE , 

.15 

46 UZC 

9.9 


.58MLF--: 

: Mo - 



47 MAM 

99 

■ 

59 SRU . 

20 '■ 




DHU5 NOTE: TT » DE05 C| RRRR, T5N7HS OF.-MMj 99 S 999-9, , MISSING 


WWEE SMUS S'-SAUS 26 OCT 75 


ID SJN XXX NN 

TT 

RRRR 

ID- STN ■ 

XXX NN 

IT 

RRRR 

67 290 

59' 

0 • 

74 270 ■- 


60 

0 

-68 IPL 

64 

0 

. ■75 MRF . ■ ' 


52 

0 

69 280 

64 

■0 

>■-76.261. • 


59 

0 

70 274 

64 

0 

77 COT 


63 

0 

-71..FHU - 

60 

0 

78 MFE- : • 


61 

0013 

72 DUG 

58 

0 

' 79 250 


60 

T 

C;73 .DMN"4.M ->i.. ' . 

56 

0- 

•- - - _ 




TEMPS DEGMbS3?'1? 

; ■ . ' ••■ • 

•1 

iiiH 

• • t r 

UNokfeolTHS OF 

IN.l. 

>- 




:.»Uf«JAY 26 ,QCf ^973/ 
^SaY 1 of WEEl^- 35 > 


iORIGINAI} PAGE IS 

€F POOR QUALITY ^ 


RF<RF 


0000 

0000 

0000 

9999 

0132 

0000 

0000 

0000 

0000 

0004 

0000 

0070 


DATA, 



temperature is reported for each meteorological station under 
the TT heading. Missing data is indicated by 99. 

The surface meteorological data for some stations are contained 
in the meteorological station statistics printout. For all 
stations in the prime verification network, the surface data 
are reported under TMET. Some care must be used when extracting 
surface data from this source. TMET will be the surface station 
data only for the ground truth control stations. 

A list of these stations is given in table 2-V. For any station 
not in this list, the only source of data is the daily National 
Weather Service report. As an example, the data for Brownsville, 
Texas for 15-29 March 1975, has been extracted and is shown in 
table 2-VI. 

2. 1.2.2 Method of Analysis 

The raw data is used to determine the appropriate thermal 
thresholds as described in section 2.3 and to calculate the 
rainfall coefficients as described in section 2.2. The results 
of those analyses are used to screen the raw radiometric data 
for cloud contamination. 

For the 15-29 March 1975, time period used in the previous 
example, the following thermal thresholds were used: 

• Day Threshold - 289. 5°K 

• Night Threshold = 280. 0°K 

Any radiometric data less than the threshold is considered 
cloudy and is deleted. 
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TABLE 2-V. - GROUND TRUTH ZONE CONTROL STATIONS 


STATION NAME 

SEDS ID 

LATITUDE 

LONGITUDE 

ELEVATION 

Hermos i 1 1 o 

01 

29:04 

110:57 

211 m 

Chihuahua 

02 

28:38 

1 06:04' 

1 354 m 

Guaymas 

03 

27:55 

110:53 

1 6 m 

La Paz 

07 

29:04 

110:20 

10 m 

Durango 

09 

24:01 

104:40 

1889 m 

Mazatl an 

10 

23:1 2 

106:25 

3 m 

Ciudad Victoria 

11 

23:44 

99:08 

321 m 

Soto la Marina 

12 

23:46 

■98:1 2 

25 m 

Tampico 

14 

22:13 

97:51 

12 m 

Topic 

1 5 

21 :30 

104:53 

915 m 

Rio Verde 

17 

21 :56 

100:00 

991 m 

Guadalajara 

18 

20:36 

103:23 

1 567 m 

Pac huca 

19 

20:08 

98:44 

2426 m 

Tuxpan 

20 

20:57 

97:23 

14 m 

Mexico City 

23 

19:24 

99:11 

2298 m 

Puebl a 

25 

19:03 

98:1 2 

2162 m 

Veracruz 

26 

19:09 

96:07 

13m 

Oaxaca 

30 

17:03 

96:43 

1 500 m 

Acapul CO 

31 

16:50 

99:55 

82 m 

Sal i ria Cruz 

32 

16:10 

95:12 

70 m 

Las Casas 

33 

16:44 

92:38 

2128 m 

Tucson 

•70 

32:07 

110:56 

778 m 

Fort Huachuca 

71 

31 ;35 

110:20 

1426 m 

Douglas 

72 

31 :27 

109:36 

1250 m 

Demi ng 

73 

32:15 

107:42 

1316 m 

El Paso 

74 

31 :48 

106:24 

1192 m 

Marfa 

75 

30:22 

104:01 

1479 m 

Del Rio 

76 

29:22 

100:55 

313 m 

Cotulla 

77 

28:27 

99:1 3 

146 m 

McAl 1 en 

78 

26:11 

98:14 

34 m 

Brownsvi 1 1 e 

79 

25:54 

97:26 

6 m 

Midi and-Odessa 

83 

31 :57 

1 02:1 1 

868 m 

San Angelo 

84 

31 :22 

100:30 

579 m 

Abi 1 ene 

85 

32:25 

99:41 

543 m 

San Antonio 

93 

29:32 

98:28 

240 m 

Austin 

94 

39:18 

97:42 

182 m 

Waco 

95 

31 :37 

97:13 

153 m 

Corpus Christi 

96 

27:46 

97:30 

12 m 

Victoria 

97 

28:51 

96:55 

32 m 1 

Houston 

98 

29:58 

95:1 2 

29 m 1 

Port Arthur 

99 

29:57 

94:01 

5 m 


2-13 





TABLE 2 -VI. - SAMPLE RADIOMETRIC DATA AS PREPARED FOR AMLYSIS 


Station Name: 

Brownsville, 

Texas 



Station ID: 

79 




Station Altitude: 6 m 




DATE 

TSDK 

VIS 

TSNK 

TT 

3-15-75 

270.5 

68 

295.5 

293.15 

3-16-75 

243,5 

82 

- 

294,27 

3-17-75 

290.5 

66 

- 

297.60 

3-18-75 

294.0 

40 

- 

298.71 

3-19-75 

311.0 

35 

295.0 

292.04 

3-20-75 

301.0 

39 

- 

293.15 

3-21-75 

302.0 

52 

- 

297.04 • 

3-22-75 

- 

- 

- 

295.93 

3-23-75 

298.5 

97 

302.0 

299.82 

3-24-75 

- 

- 

- 

298.15 

3-25-75 

283.0 

34- 

- 

292.60 

3-26-75 

251 .0 

99 

- 

298.71 

3-27-75 

238.5 

74 

234.0 

299.82 

3-28-75 

- 

- 

- 

300.93 

3-29-75 

269.0 

113 

254.5 

285.38 
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The rainfall coefficient defines a decision field which assigns 
all possible combinations of TSD and VIS to one of five classes. 
Only class 1 is cloud free. A part of the decision field used 
during the time period of the example. is shown in figure 2-3. 

Consider the data for 17 March 1975: 

TSDK = 290.5 
VIS. = 66 

The value of TSDK is greater than the day threshold of 289.5 but 
the VIS, TSDK combination entered in figure 2-3 does not fall 
into the clear category. Table 2-VII shows the usable data that 
remained after thresholding and cloud discrimination- In this 
way, a data set can be assembled suitable for regression analysis. 

The initial data sets should contain at least 50 points for each 
of the three cases. If insufficient data are available for one 
or more cases, either the time period or the number of stations 
should be expanded. 

Multiple linear regression analysis assumes a model of the 
following form: 

y - 3 + + ••• + 3 y 

o 112 2 n n 

The data are collected and may be written as 

^17 •’ ^21 ^ ' ’ ' * ^ 1 

^12 ^ ^22 ^ ‘ ^ ^ n 2 ^ ^2 


Im 2m ^ nm m 


7.-T 5 



TSDK 
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Figure 2-3. — Cloud discrimination decision field for March 1975 data 


to M ^ J I\3 N5 N>( O ro N) to 1 J I O to r\3 hJ fO N3 to to to to Ivj to IvJ ro ts) to ^0 hJ M to fO t'J to O t J f- J to h J h. 


TABLE 2-VII. - SAMPLE RADIOMETRIC DATA FROM TABLE 2-VI 
AFTER THERMAL THRESHOLDING AND CLOUD DISCRIMINATION 


Station 

Name : 

Brownsville, Texas 

Station 

ID: 

79 


Station 

Altitude : 

6 m 


DATE 

TSDK 

TSNK 

TT 

3-15-75 

— 

295.5 

293.15 

3-16-75 

- 

— 

_ 

3-17-75 

— 

— 

— 

3-18-75 

— 

— 

— 

3-19-75 

311.0 

295.0 

292.04 ' 

3-20-75 

301.0 

— 

293.15 

3-21-75 

— 

— 

— 

3-22-75 

— 

— 

— 

3-23-75 

— 

302.0 

299.82 

3-24-75 

— 

— 

— 

3-25-75 

— 

— 

— 

3-26-75 

— 

— 

— 

3-27-75 

— 

— 

— 

3-28-75 

— 

— 

— 

3-29-75 

— 

— 

— 
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Assume that the I. are measured exactly, but that 7 has 
random, normally distributed error with constant variance. 

Least squares analysis provides a method of finding estimates 
h . of the parameters 0. . These estimates are the regression 

'V 'ly 

coefficients . 


The procedure is to minimize the sum of squares 
S = Z (Y - 0 


0 


by setting ^ = 0. ^ 


0 


BfXi 02X2 - .... - 


= 0, and replacing 0 q, 0^ by their 


estimates bg,.b^ . 


Unsubscripted summations are over the number of observation 


II 

.H 

2 , , . . n and j = 1 , 

2 , . . . , 

n. 

Defining, 

c . . = Ex . X - 
ij 1 1 

= EX.X. 
1 D 

EX. EX. 

1 i 

m 




EX. EY 


^iy = ^^i^ 

= EX^Y 

1 

m 


c = Ey^ 

yy 

- 

CEY)^ 

m 

leads to 

the n simultaneous 

linear 

normal equations 


Ciibi + <^12^2 * 

• • • ^ ^In^n ^ly 


^ 21^1 ‘^ 22^2 


+ c„ b 
2n n 


= c 


2y 


I 


c lb. + c „b,, + 
nl 1 n2 2 


+ c b 
nn n 


= c 


ny 
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Together with 


bo = Y-Zb.X 
1 


The system I is conveniently solved by inverting the n x n 

coefficient matrix C = (c..) to give the matrix C ^ = Cc"-) 

1 J 1 J 

(The program checks to see i£ C is singular and if so, it 

will print out.) We then have 


b . = Z cl . c . 

1 j iJ jy 


The variance o is estimated by 


2 

s 


(c - Z b • c . ) 
^ yy ^ ± ly^ 

(m - n - 1)' 


there now being only (m - n - 1) degrees of freedom since 
(n + 1) coefficients have been estimated from the data. The 
variances and covariances of the b^ are given by 

var (b^) = s^ c^^ 
cov Cbj_,b^^) = s^ cr^ . 

It may also be shown that 

var (b^) = s^/m + s^ Z Z X. X. cl. 

0 . . 1 j 11 

1 J J J 

and 


cov CbQ,b^) 


2 — 
s Z X. 


c . . 

ij 
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is computed as 


2 

The coefficient of multiple determination R 
folloxsTS : 


2 f i ‘'iy 

yy 


Its significance, and hence the overall significance of the 
regression, is tested by analysis of variance, and thus by 
obtaining the ratio 


_ , (m - n - 1}R^ 
n(l - 

and testing it against the F-distribution with n and 
m - n " 1 degrees of freedom. 

An adjusted coefficient of multiple determination, corrected 
for sample size and number of independent variable often useful 
in comparing results, is computed as follows: 

r2 = 1 _ (m - 1)(1 - 
m - n - 1 

The quality of the final regression fit can be improved by 
examining the residuals from the first run for each case. In 
general, any data set which has an error greater than twice 
the standard error of estimate should be deleted and the 
regression rerun. The final coefficient set can be obtained 
in, at most, three passes. 

The correlation matrix, multiple correlation coefficients, 
standard error of estimate, analysis of variance, and regression 
coefficients should be retained for future reference. 
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Table 2-VIII shows the results of the final regression analysis 
for the 15-29 March 1975 SEDS night only radiometric data. In 
this example, the resulting regression is highly significant 
(F-Value = 179.1) with an excellent standard error of estimate of 
1.90°C. The correlation matrix shows strong correlation between 
mean air temperature and both of the independent variables. When 
evaluated over all data during the 15-29 March 1975 time period, 
the performance of this coefficient was very good with a root 
mean square (rms) error of 2.41°C. 

The SEDS DMAT coefficients are found by transforming the 
regression coefficients from equations (2-1) to (2-3) into the 
SEDS computational form shown in equations (2-4) to (2-6) . The 
transformations are given in equations (2-7) to (2-16). Table 
2-IX gives the coded and uncoded coefficients as found in the 
15-29 March 1975 DMAT calculation. 

2.1.3 VERIFICATION 

The verification of the coefficients ’ set is a continuing 
process ultimately resulting in the recalculation of the DMAT 
coefficients. The performance of the DMAT estimates should be 
evaluated over the prime verification network shown in 
table 2-III. 

The verification uses a measure of the mean and standard devia- 
tions of the observed errors. The error is defined as 

ERROR = Observed - Predicted 

The Bias error is then the arithmetic mean of the errors. A 
positive error means that SEDS underpredicted with the DMAT 
estimate. The rms error is the standard deviation of the errors. 
Approximately 67 percent of the values will be within one 
standard deviation of the mean and 95 percent will be within 
two standard deviations. 
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TABLE 2-VIII. - MULTIPLE REGRESSION RESULTS FOR 15-29 MARCH 1975 
FOR CASE 3 - ONLY NIGHT RADIOMETRIC DATA USABLE 


Correlation Matrix: 



TSNK 

ALT 

TT 

TSNK 

1.000 

1 

c^ 

o 

.968 

ALT 

- .670 

1.000 

- . 714 

TT 

.968 

- .714 

1.000 


Multiple Correlation Coefficient : .972 

Standard Error of Estimate : 1.901°C 

Analysis of Variance for the Multiple Linear Regression : 



Degrees 

Freedom 

Sum 

Squares 

Mean 

Squares 

F 

Value 

Due to Regression 

2 

1294.86217 

647.43108 

179.12646* 

Deviation about Regression 

21 

75.90197 

3.61437 


Total 

23 

1370.76415 



*Signif leant at the .99 level. 




Regression Coefficients: 
Constant = 47.89565 
X(TSNK) = .83225 

X(ALT) = -.00155 







TABLE 2-IX. - 15-29 MARCH 1975 DMAT COEFFICIENTS 
IN BOTH CODED AND UNCODED FORMS 


Coded ; 

DMAT = AOC + 
DMAT = BOC + 
DMAT = COC + 

AOC = -126.69744 
AIG = 001.23440 
A2C = -000.13490 
A3C = 000.01000 

Uncoded : 

DMATK - 
DMATK = 
DMATK = 

= 90.03634 
o 

A^ = 00.61720 
Ag = 00.06745 
A^ = 00.00010 


AlC (TSN) 
BIC (TSD) 
CIC (TSN) 

+ A2C (TSD) 
+ B2C (ALT) 
+ C2C (ALT) 

+ A3C 

(ALT) 


BOC = 

-82.66244 

COC = 

-175. 

0970 

BIC = 

01.07400 

CIC = 

001. 

6651 

B2C = 

-00.19000 

C2C = 

-000. 

1550 


A^ + A^TSNK + AgTSDK + A^ALT 

B + B.TSDK + B„ALT 
O 1 2 

C^ + C-TSNK + C„ALT 
o 1 2 


B = 130.67039 

o 

B^ = 000.53700 

= -000.00190 


C = 47.89565 

o 

C^ = 00.83225 

C^ = -00..00155 
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The Bias and rms errors for the original data set upon which the 
DMAT coefficients were calculated should be found- These values 
form the base line to which future errors should be compared. 

In the example from 15-29 March 1975, the Bias error for Case 3 
(night radiometric data only) was +.52®C and the rms error was 
2.41°C, By 10-23 May 1975, the Bias error had changed to 1.43°C 
and the rms error to 3.63°C. At that time the DMAT coefficients 
were recalculated. For the same case, the new Bias error was 
0.0°C with a rms error of 2.16°C. Thus, providing the base-line 
figures for the next evaluation. 

In general, the errors for at least a week's data should be 
calculated to avoid random sampling errors. When either of the 
errors has increased more than 1.0°C, the DMAT coefficients 
should be recalculated. 

2.2 RAINFALL COEFFICIENTS 

The primary method used by SEDS for detecting cloud contaminated 
day radiometric data is the rainfall classifier. The day pass 
infrared and visible data are classified into one of five classes 
by a set of discriminant functions. Only the data which falls 
into the "clear" category is used by SEDS for calculating the 
daily mean air temperature estimates. 

2.2.1 BACKGROUND MATERIAL 

The name of the rainfall classification algorithm is misleading. 
As used by SEDS it is actually a cloud classification algorithm. 
This feature of SEDS was intended to serve the dual purpose of 
providing information about cloud contamination in the radio - 
metric data and of providing a research tool for making estimates 
of precipitation. The research into estimating precipitation 
probability or amount from SEDS data .was not carried out. The 
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algorithm provides a reliable means o£ eliminating cloud contami- 
nation from the day radiometric data. 


2 , 2 . 1 . 1 Theory 

Briefly stated, the problem is one of differentiating between two 
or more populations on the basis of multivariate measurements. 
Given the existence of N populations, there exists a two-part 
problem. The first is to correctly distribute a set of samples 
among the N populations; the second is to determine a function 
based on these samples, which will allow the assignment of a new 
sample to the correct population. The first class of problems 
is generally called classification analysis, and the second is 
called discriminant analysis. 

Classification analysis has been highly developed for applica- 
tions as diverse as numerical taxonomy and land use classifica- 
tion. The problem is to classify a series of samples into groups 
which are as distinct as possible. 

In SEDS , the day pass data consists of a pair of multispectral 
measurements for each point. The infrared and visible data when 
used as a pair provide more information on cloud' cover than 
either data could by itself. Through classification analysis, 
it is possible to distribute a series of samples into the most 
appropriate groups. By controlling classification procedure, it 
is possible to create a single cloud-free class and several 
cloudy classes. 

Discriminant analysis takes the samples placed in each group and 
derives a set of functions that enables the assignment of any 
new sample to the group which has the highest probability of con- 
taining the sample. These functions serve as the cloud classifi- 
cation algorithm in SEDS. The resulting decision field has been 
shown previously in figure 2-3. 
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2. 2. 1.2 SEDS Utilization 


The algorithm used in SEDS is a maximum likelihood classifica- 
tion using linear discriminant functions. The algorithm assigns 
any possible input sample to one of five class regions or to an 
"unclassif iable" class. A thresholding capability allows the 
operator to set the probability level for inclusion in a class. 
To be used in the DMAT calculation an input sample must be 
classified clear. The classification algorithm may be defined 
as follows: 


II 

^10 

+ 

a^^VIS 

+ 


(2-17) 

^2 = 

^20 

+ 

a^^VIS 

+ 


(2-18) 

Ks - 

^30 

+ 

a3^VIS 

+ 

ajzTSD 

(2-19) 

II 

^40 


a^^VIS 

+ 

a^^TSB 

(2-20) 

^5 ” 

^50 

+ 

a^^VIS 

+ 

a TSD 

(2-21) 


Where the symbols are defined as follows: 

a.. - The rainfall estimation coefficients 

VIS - Eight-bit value of the visible detector voltage 

TSD - Eight-bit coded value of the day radiometric temperature 

The sample is assigned to the class which has the largest 
for the (VIS, TSN) vector of that sample. Through the use of a 
threshold value, it is possible to emphatically set the probabil- 
ity level for inclusion in each class. When a given vector's 
probability of inclusion in any class is low, the vector is 
assigned to a null class. To date, this feature has not been used. 
Although the classes have no meaning beyond the clear/cloudy 
classification, the full results are given in the Meteorological 
Station Statistics Printout under "R CODE”. Table 2-X gives the 
interpretation of each of the symbols in the printout. 
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TABLE 2-X. - METEOROLOGICAL STATION STATISTICS 
PRINTOUT R CODE INTERPRETATION 


Symbol 

Class 

Interpretation 

C 

Dry 

Clear 

CD 

Cloudy Dry 

Cloudy 

L 

Light Rain 

Cloudy 

M 

Moderate Rain 

Cloudy 

H 

Heavy Rain 

Cloudy 

0 

Null Class 

Cloudy 


For computational purposes, SEDS requires that the discriminant 
coefficients a^^ be scaled. The scaled coefficients, , are 
input to SEDS via the SEDS System Parameter Control Form shown 
in figure 2-1. The coefficients are scaled integers subject to 
the following constraints: 

• The largest multiplicative term £ 255 

15 

• The absolute value of the largest constant term _< 2 
2.2.2 METHOD OF CALCULATION 

The approach consists of collecting a set of data from SEDS which 
is classified into five different groups by the classification 
program ISOCLS/PRINT . It identifies groups or clusters in the 
input data that have similar elements within clusters that are 
different from the elements of other clusters. The program uses 
builtin default values for all classification parameters. The 
samples from each group are then used to calculate the coeffici- 
ents for equations 2-17 through 2-21. The discriminant function 
coefficients are calculated by the program DISCRM, 

Program tapes and listings of these programs along with user docu- 
mentation are available from Dr. C, M. Barnes, Health Applications 
Officer, Life Sciences Directorate, NASA, JSC. Details of these 
programs will not be repeated here. 
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2 . 2 . 2 . 1 Data Used 


The radiometric data source is the SEDS meteorological station 
statistics printout. A typical printout was shown previously in 
table 2- II. The paired values for TSDK and day VIS for all 
stations are extracted. In general, two weeks of SEDS data are 
required to be sure of obtaining a sufficient amount of data. 

The data are key punched for the program ISOCLS/PRINT in the 
following format: 

a Col 1-3 — Sequence number of the data card 1, 2, 31.. (not 
required by program) 

a Col 9-12 — TSDK as four digits without decimal 
a Col 25-27 —VIS as three digits without decimal. 

The primary source of data for verification is a report of sky 
conditions supplied by the National Weather Service. 

2 . 2 . 2 .2 Analysis Used 

The analysis used is a two-step procedure. The raw data is 
classified into five classes by ISOCLS/PRINT. The samples in 
each class are then analyzed by DISCRM to provide the unsealed 
discriminant coefficients. 

The program ISOCLS/PRINT performs an unsupervised classification, 
often called clustering, of the input data. The program will 
accept as input a set of cluster means as initialization values. 
The number of iterations required is greatly reduced by this 
initialization, and in this case the quality of the final result 
is better . The input cluster means are one of the results of the 
previous calculation of classification coefficients. 

For example, the analysis for the spring evaluation was . performed 
on 15-29 March 1975 SEDS data. The program ISOCLS/PRINT was 
initialized with the following set of class means which had to 
be calculated for late summer data. 
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MEANS 


Cluster 

Visible 

Infrared 

1 

36.34 

304.44 

2 

44.52 

291.26 

3 

70.42 

264.58 

4 

115.60 

278.77 

5 

116.17 

256.83 


Table 2-XI shows the statistical results o£ ISOCLS/PRINT after 
six iterations. As can be seen, the class numbers have been 
reordered and the class means have in some cases undergone con- 
siderable change. The clear class is the one with the highest 
TSDK value and the lowest VIS value. As TSDK decreases and VIS 
increases, the classes represent successively higher and more 
dense clouds. 

ISOCLS/PRINT also outputs a list of the members of each class. 

The data cards are reordered on the basis of this list and input 
to the program DISCRM. DISCRM corrects the TSDK data to the 
eight-bit coded TSD used in SEDS. The VIS data is unchanged. A 
sample of the program output is shown in table 2-XII. The mean 
scores are the resulting class means for each variable. The dis- 
criminant coefficients are given for each class. Additionally, 
the program produces a table of the classification results for 
all input data. If the ISOCLS/PRINT program worked correctly, 
only a few outlying points will be classified incorrrectly . 

The unsealed and scaled coefficients as used by SEDS are shown 
in table 2-XIII. The scale factor is found by dividing the 
largest multiplicative term into 255. In this example 255/. 98661 
= 258.46079. This algorithm is the primary cloud detector for 
the day radiometric data. A portion of the decision field has 
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TABLE 2-XI. - TYPICAL STATISTICS PRODUCED BY 
ISOCLS/PRINT FOR 15-29 MARCH 1975 SEDS DATA 
(SELECTED STATIONS ONLY) 


intermediate Printout for iteration 


TOTAL number Op CLUSTERS 3 B 
total number Of points = 10s 


CLUSTER 

1 

2 

3 

5 


symbol 

POINTS IN CLUSTER 

1 

26 

2 

9 

3 

18 

8 

1 3 

5 

39 


means 


cluster 

1 

2 
3 
H 
5 


CH( !) CH{ 2) 
296,25 S7.DS 

261,28 59,22 

278,75 118,56 

235,77 107.85 

303.22 38.78 


standard Deviations 


cluster 

1 

2 

3 
8 
5 


CH ( t ) 

CHI 2) 

6 » 00 

7,95 

12.96 

7,55 

12.00 

15.98 

1 0. 59 

20 . 2 1 

6.43 

8.3** 


cluster 

1 

2 

3 

8 

5- 


DIsTANcES -between clusters 


1 

2 

3 

4 

5 

.00 

4,08 

6.27 

8 .59 

2.35 

S.D8 

,00 

5.22 

4,50 

6 . 27 

6,27 

5,22 

tOG 

3,7 1 

9 . 95 

8,59 

4,50 

3.7 1 

• 00 

11.30 

2,35 

6,27 

9,95 

11.30 

. 00 


Driginac page is 

OF POOR QUALITY 
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TABLE 2-XII. - SAMPLE OUTPUT OF DISCRJI FOR lS-29 MARCH 1975 SEDS DATA . 




GROUP 

1 

2 

3 



SAMPLE. 

2.6 

9 

1 rt 


1 3 

MEAN SCORES 
1 

18a, 5COOG 

1 1 B • 55556 

lS3,50DOO 

67 

.538^6 

2 

'17,07692 

59,27222 

t IH, 55556 

107 

,8'16 15 


39 

202.i»359c 

3<*.7H359 


Total 
1 05 


FUNCTION 

1 

2 

3 

1 

5 

coefficient 

I 

.57092 

.37160 

.50186 

.25223 

. 6Q165 

2 

, R 7 1 17 

•53359 

.96661 

,68819 

• 3809 t 

constant 

-6H.90C30 

-37.91660 

-95.25861 

-56,12763 

-67.83019 


CLASSIF IC A r i ON 

MaTX 1 X 


FUNCT I ON 

1 

2 

GROUP 



1 

23 

n 

2 

1 

8 

3 

T 

C 

1 

rs 

1 

5 

1 

1 


3 

1 

s 

TOTAL 

0 

c 

3 

26 

0 

0 

r< 

9 

18 

0 

n 

16 

-3 

1 2 

0 

13 

0 

n 

37 

39 



already been shown in figure 2-3. The uses of this field are 
also discussed in section 2. 1.2.1. 


TABLE 2 -XI II. - SCALED AND UNSCALED SEDS 
RAINFALL CLASSIFICATION COEFFICIENTS 


Unsealed 


O 

II 

-67.83849 

^30 " 

-37.91668 

^50 

-95.25861 

^11 = 

.38091 

^31 " 

, 53059 

^51 = 

. 98661 

^12 

. 60485 

^32 

.37460 

^52 

. 50486 

^20 " 

64.90030 

^40 " 

-56.42783 



^21 

.47117 

^41 

,88849 



^22 

. 57092 

^42 

.25223 



Scale 

factor : 

258.46079 






Scaled 



> 

1-* 

o 

II 

-17533 

^30 

-98000 

^50 

-24621 

^1 = 

98 

^31 

137' 

^1 = 

255 

^12 " 

157 

^32 

97 

^52 " 

130 

^20 

-16774 

^40 

-14584 



^21 

122 

^41 

230 



^22 

148 

^42 " 

65 
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2.2.3 VERIFICATION 

SEDS produces an image of the results of the rainfall classifi- 
cation algorithm. This product is extremely useful in verify- 
ing the quality of the classification. CTear and cloudy areas 
can often be identified with certainty in the isothermal imagery. 
Any conflict with the rainfall classification is apparent. 

One of the key problem areas with this algorithm has been its • 
classification of the highly reflective desert region. Often an 
otherwise satisfactory set of coefficients will consistently 
classify the bare soil and sand of the desert as cloud. The 
rainfall image is vital to identifying this problem when it 
occurs . 

Numerical verification depends upon surface observations of cloud 
cover. This type of observation is difficult to interpret since 
it is an integrated average of what the observer sees over the 
entire hemisphere of the sky. On the other hand, the radiometer 
views a single kilometer square and may well find the only cloud 
in an otherwise clear sky. 

Table- 2-XIV shows an example of the verification of performance 
of this algorithm with surface data. As shown, the performance 
of the algorithm can be considered satisfactory. This is an 
example j of course, any real evaluation would require more data. 


2-33 



TABLE 2 -XIV. - COMPARISON OF SURFACE REPORTS 
FROM TEXAS WITH SEDS CLOUD DETERMINATION 
FOR 5-11 JUNE 1975 DATA 


Surface cloud cover 
reported in lOths 

SEDS Classification 
Clear Cloudy 

0-2 

7 

O’ 

3-5 

3 

0 

6-8 

2 

7 

9-10 

2 

27 


2.3 THERMAL THRESHOLDS 

The SEDS thermal thresholds are . temperature limits used to screen 
radiometric data for cloud contamination. Any radiometric tem- 
perature less than the applicable thermal threshold is not used 
in the calculation of DMAT estimates as described in section 2.1. 
The thermal thresholds for the day and night radiometric data 
are independently calculated and used in SEDS. 

2.3.1 BACKGROUND MATERIAL 

The problem is to identify cloud-free radiometric data on the 
basis of the temperature of the emitting surface. This approach 
makes use of the fact that the temperature of the lower atmos- 
phere decreases with increasing height above the earth's surface. 
Since the tops of clouds are at the temperature of the air at 
that level, it is reasonable to expect the radiometric tempera- 
ture to be less when viewing clouds than when viewing the ground. 

The infrared radiometric data processed by SEDS shows tempera- 
tures ranging from less than 230°K to over 320'’K. Figure 2-4 
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FREQUENCY 



Figure 2-4 




shows a histogram o£ the day radiometric temperatures over 
meteorological stations in the SEDS image during the 1-10 August 
1975 time period. The coldest temperatures are from the tops of 
the highest clouds, while the warmest come from the desert 
regions of Mexico and the southwestern U.S. 

However, the differences between cloud top and surface tempera- 
tures are not always great. In addition, the surface tempera-^ 
ture in mountainous areas may well be less than the temperature 
of clouds occurring at lower altitudes. Careful analysis of 
histograms of the radiometric data provide estimates of the 
thermal thresholds. 

SEDS software compares each pixel of data with predetermined 
temperature limits. Any radiometric temperature value greater 
than the threshold value is considered to be cloud free. The 
day thermal threshold is generally several degrees warmer than 
the night threshold. The threshold values change seasonally 
with the highest values occurring during the fall. All changes 
in the thermal thresholds are initiated via the SEDS System 
Parameter Control Form shown previously in figure -2-1. 

2.3.2 METHOD OF CALCULATION 

The basic method of calculation uses a histogram constructed 
with data taken from the SEDS meteorological station statistics 
printout. Fine tuning of the threshold values may be carried 
out by examining the residual errors of the DMAT estimates. 

A given set of thermal thresholds remain valid for 4-6 weeks. 
Recalculation is generally carried out in conjunction with 
determination of the DMAT coefficients. 
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2. 3. 2.1 Data Used 


The data source is the daily SEDS meteorological station statis- 
tics printout shown previously in table 2 -II. The TSD and TSN 
data are extracted for all stations over a one-week period. 
During extremely cloudy periods, it may be necessary to use two 
weeks of data in order to have a significant number of cloud- 
free points. 

2.3. 2. 2 Analysis Used 

An analysis of histograms of the raw radiometric data is the 
approach used. Complete histograms are not required since 
extremely cold temperatures may be assumed to represent clouds. 
Separate histograms are plotted for the day and night data. 
Figures 2-5 and 2-6 show examples of these plots for the 
28 September to 9 October 1975 time period. 

The best thermal threshold value has been found to' be 5-7 
degrees less than the first major frequency maximum occurring 
for temperatures greater than 280“K. Table 2-XV shows these 
values for our example. 

Table 2-XV. - TEMPERATURE OF FIRST CLOUD- FREE FREQUENCY MAXIMUM 
AND ESTIMATED THERMAL THRESHOLD FOR 28 SEPTEMBER 
TO 9 OCTOBER 1975 SEDS DATA. 



Temperature of 

Thermal 


Frequency Maximum 

Threshold 

Day 

303° 

296° 

Night 

291° 

285° 


The histogram itself indicates the correct choice by a noticeable 
discontinuity. 
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FREQUENCY 



Figure 2-5. — Histogram of day radiometric data 
28 September - 9 Octoljer 1975. 




Figure 2-6 — Histogram of night radiometric data 
28 September - 9 October 1975. 
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2.3.3 VERIFICATION 


The evaluation of the thermal thresholds are linked to the 
evaluation of the DMAT estimates as described in section 2.1.3. 
When large rms errors are noted in the DMAT estimates, the 
radiometric data should be histogrammed in order to determine if 
the thermal thresholds should be changed. 

When the temperatures are increasing in the spring through late 
summer, an analysis o'f the residual errors associated with" the 
DMAT estimates provides additional information. When the 
residual error is consistently large for temperatures just above 
the existing threshold, it is a clear indication that recalcula- 
tion is necessary. 

2.4 COMBINED SCREWWORM PRODUCT COEFFICIENT 

The SEDS combined screwworm product is a relatively late develop- 
mental product. Its purpose is to provide a single test estimate 
of the potential growth of a generation of screwworms as a func- 
tion of weather conditions. SEDS was designed to provide four 
separate estimates of the potential population growth detailing 
the effect of temperature- on adult and pupae survival, of soil 
moisture on pupae survival, and of mean pupation time. The 
combined product was implemented to provide a regression based- 
combination of the four potential growth estimates. 

2.4.1 BACKGROUND MATERIAL 

This section is -intended to give an overview of the development 
leading up to the combined screwworm product, which is still 
very much in the developmental phase. Although the mathematical 
form of 'the algorithm implemented in SEDS has been statis- 
tically verified, the regression coefficients are derived from 
simulations based on historical surface data. In addition, the 
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well !known limitations on applying regression results outside o£ 
their developmental area make the general applicability o£ the 
combined product a matter for further research. 

2.4. 1.1 Independent Variables 

There are four variables calculated by SEDS that serve as the 
basic information from which the combined screwworm product is 
calculated. These variables are recalculated and updated daily 
for each pixel in the 625 element by 550 line data base array. 
The following is a brief summary of the manner in which SEDS 
handles each variable. 

The short term mean air temperature (STMAT] depicts the effect 
of daily mean air temperature on the growth and survival of the 
adult screwworm fly. Computationally, the STMAT is determined 
by the following low pass filter. 

STMAT = .44 DMAT + .56 [STMAT (-1)] 

where DMAT - Eight-bit coded values of the 

daily mean air temperature 
estimate as described in 
section 2 . 1 . 1 . 2 . 

STMAT (-1) " The STMAT value calculated 
the previous day. 

The long term mean air temperature (LTMAT) shows the effect of 
daily mean air temperature on the growth and survival of the 
screwworm pupae. Computationally, the LTMAT is the running sum 
of the DMAT estimates over the last N days-. 

The crop moisture index (CMI) is a measure of soil moisture 
deficit as described in section 3.1. The CMI variable used in 
SEDS is an N-day sum of the daily CMI values. This variable is 
intended to depict the effect of soil moisture on the screwworm 
pupae . 
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The degree day sum (DDSUM) variable is an estimate of the screw- 
worm's mean pupation time. This variable is used directly as an 
independent variable and to estimate generation length. Com- 
putationally, the DDSUM is an N-day accumulation of the screw- 
worm degree day values. 

N 

DDSUM = I (DMATCi) “ 284] when DMAT(i) > 284°K. 

L=1 

The application of the degree day sum to screwworms is given 
elsewhere (ref. 1). 

The STMAT, LTMAT, CMI , and DDSUM values are subject to additional 
calculation and coding prior to the determination of the combined 
product . 

The STMAT is converted to a growth potential value using a look- 
up table based on figure 2-7. The development of this figure is 
given by Phinney (ref.l] and discussed briefly in section 1.1. 

The growth potential is encoded by 

= integer [(100 In GP) + 128.5] 

\ 

The LTMAT sum is divided by N to determine the mean and con- 
verted to a growth potential again using figure 2-7. The growth 
potential is encoded by 

X2 = integer [(100 In. GP) + 128,5] 

The CMI sum is divided by N to determine the mean and con- 
verted to, a growth potential using a look-up table based on 
figure 2-8. The growth potential is encoded by 

X^ == integer [(100 In GP) + 128.5] 

The DDSUM is used to calculate mean pupation time (MPT) using 
the following -formula. 

MPT = NC139.1/DDSUM) 
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GROWTH POTENTIAL 



Figure 2-7.— Relation between mean air 
temperature and growth per generation. 


GROWTH POTENTIAL 



MEAN SOIL MOISTURE AS MEASURED BY 
THE PALMER CROP MOISTURE INDEX 


Figure 2-8.-= Relationship between CMI and growth potential. 



The estimated generation length is MPT + 15. The generation 
length is converted to a growth potential value using a look-up 
table based on figure 2-9. The growth potential is encoded by 

= integer [(100 In GP) + 128.5] 

The variables available for calculating a combined screwworm 
growth potential are 

- Eight-bit coded growth potential 
calculated from the STMAT 

X 2 - Eight-bit coded growth potential 
calculated from the LTMAT 

■ X^ - Eight-bit coded growth potential 

calculated from the CMI 

X^ - Eight-bit coded growth potential 
calculated from the DDSUM 

2. 4. 1.2 Model Selection 

The selection of a multiplicative model for the SEDS combined 
screwworm model was based on computation and statistical expedi- 
ency as well as mathematical plausibility. Statistically, it 
was desirable to have an intrinsically linear model so that 
multiple linear regression analysis could be used. The multi- 
plicative model may be written in the following form 

a, a, 

Y = a^X^ X 3 *^X^ (2-22) 

where ct^, a^, ^ 2 ^ and are parameters and e is the 

random error. Taking the natural logarithms of Eq. (2-22) 
converts the model into a linear form. 

5,nY = 5,nct + a, £ X- + a„£nX„ + a-£nX_ + a.£nX. + £ne (2-23) 

0 Inl 22 33 44 

The transformed model given in Eq.(2-23) is linear. Linear mul- 
tiple regression may be used to calculate estimates of the 
parameters . 
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Figure 2-9.- Relationship between generation length 

and growth potential. 




Computationally, this form was desirable in that the data base 
values to be used as independent variables were contained in 
SEDS in logarithmic form. 

2,4.1. 3 Model Verification 

In an effort to determine the validity of the mathematical form 
of the combined screwworm product algorithm, three basic models 
were examined. First was- the linear model which may be written 

-1- a^GPCGEN) + e 

where GP(T^) - Growth potential of the mean air 

temperature 

GP CTp). “ Growth potential of the mean air 
temperature during the pupae 
phase 

GP(CMIp) - Growth potential of the mean CMI 
value during the pupae phase 

GP CGEN3 - Growth potential of the estimated 
generation length. 

Second was an exponential model which may be written 

Y = exp [a^ + a^GP(T^) + GP(T^) 

+ a^GPCMp) + a^GPCGEN3]c. 

In linear form, this model may be written 

£nY = a + a, GP(T.) + a«GP(T ) 
o 1 A 2 p 

+ a^GPCCMIp) + a^CGEN) Me. 
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The third was the multiplicative model implemented in SEDS which 
was given in section 2.4.1. 2. 

These models were evaluated using the set of historical data 
used to develop the growth potential functions as described in 
Phinney (ref.l). Case reports less than ten per week were 
eliminated to reduce the effect of sampling error. It should be 
emphasized that although similar, the variables are not e:Kactly 
the same as the SEDS data base variables. The significance of 
these models as evaluated by the F test described in section 
2. 1.2, 2 may be seen in table 2-XVI. 

TABLE 2-XVI. - F TEST OF SIGNIFICANCE FOR THREE 
POSSIBLE SCREWWORM REGRESSION MODELS 


Model 


F 

Linear 

(4,185) 

8.5 

Exponential 

(4,185) 

15.0 

Multiplicative 

(4,185) 

16.7 


These results suggest that for relatively simple models 
utilizing these variables that the choice of SEDS- algorithm was 
correct . 

2 . 4 . 1 . 4' Developmental Phase Coefficients 

During the developmental phase of this algorithm insufficient 
SEDS data were available to derive estimates of the parameters 
for Eq.(2-23). Rather than use arbitrary coefficients, ifwas 
decided to use those derived during the model evaluation 
described in section 2. 4. 1.3. Table 2-XVII gives the coef- 
ficients derived during that study. 
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TABLE 2-XVII. - COEFFICIENTS DERIVED FOR 
MULTIPLICATIVE MODEL 


Parameter 

Coefficient 

Computed t-test value 

£na^ 

-1.58553 

• 

“l 

,37046 

2.644* 

“2 

-.17782 

-.893 

“3 

. 29818 

2.660* 

“4 

,32869 

2.050 

*Signif leant 

at the .99 level. 



For implementation into SEDS was set equal to one. Addi- 

tionally, the coefficients in table 2-XVII were scaled to take 
advantage of the full range of available colors in the output 
image. The scaling was carried out by dividing the individual 
coefficient by the sum of all four coefficients. Table 2-XVIII 
shows the coefficients used during the developmental phase. 

'TABLE 2-XVIII. -DEVELOPMENTAL PHASE COMBINED SCRE^WORM 

PRODUCT COEFFICIENTS 


Variable 

Parameter 

Coefficient 

Intercept 

«o 

1 

STMAT 

“l 

.46 

LTMAT 

“2 

-.22 

CMI 

“3 

. 36 

DDSUM 

“4 

.40 


2-49 





2,4.2 COEFFICIENT CALCULATION AND VERIFICATION 

Both, the calculation and verification of the SEDS combined 
screiworm product coefficients are currently under study. Data 
is currently being summarized iri Texas by -county for both SEDS 
and USDA case reports. When sufficient data has been accumu- 
lated for statistical significance, multiple regression studies 
will be carried out. 
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3. CALCULATION OF THE CMI 


The crop moisture index (CMI) is a measure o£ short-term anom- 
alies of soil moisture. Based on a relatively simple two-layer 
soil moisture model using weekly rainfall and temperature, the 
CMI is included in SEDS to provide a measure of the effect of 
soil moisture on the screwworm's pupae phase. 

3.1 BACKGROUND MATERIAL 

The CMI values used by SEDS are calculated by the same basic 
program used operationally in the United States by the National 
Weather Service and published in the Weekly Weather and Crop 
Bulletin . Historical data for 36 stations were analyzed to 
provide information on the normal conditions required for weekly 
operation of the CMI for SEDS. The existing program is imple- 
mented on the JSC Univac 1110 system. Operating instructions 
for the current program along with sufficient documentation for 
implementation on other machines may be obtained from C. M. 
Barnes, Health Applications Officer, Life Sciences Directorate, 
NASA, JSC. 

3.1.1 SOIL MOISTURE MODEL 

The basic CMI model assumes a hydrological balance for the soil. 
The amount of water in the soil can be calculated from a balance 
which includes precipitation, evapotranspiration, storage, 
runoff, and percolation through the soil. 

The model assumes that the first inch of water is in the upper- 
most layer of the soil, the depth of which varies from one to 
eight inches depending upon the soil type. The remaining water 
is in a second layer. Precipitation fills the first layer 
followed by the second. Amounts over the holding capacity of 
the soil are disposed of by runoff and percolation. 
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Water loss through evapotranspirat ion is calculated by the 
Thomthwaite method in "which potential evapotranspiration is a 
function of mean air temperature and day length. Losses are at 
full potential rate until the first inch of water is depleted. 
The second layer gi"ves up water at a slower rate which is an 
empirical function of the ratio of water held to the* water 
capacity. 

The Palmer Drought Index Cref. 2) was developed to give a 
measure of long-term moisture anomalies which would be useful 
for determining the scope, severity, and frequency of prolonged 
droughts. This program determines the climatological norms from 
which deviations are measured. These norms are also required 
for the CMI analysis. 

The historical analysis was performed for 32 Mexican and 4 
United States stations shown in table 3-1. The required data 
included a historical set of monthly mean temperature and total 
precipitation along with station latitude and water holding 
capacity of the soil. The set of norms can be extrapolated in a 
linear manner to nearby stations which lack sufficient climato- 
logical records to allow direct computation. 

A master data file was established on the JSC Univac 1110 for 
the stations shown in figure 3-1. This file contains the 
historical averages against which current weather is compared. 

3.1.2 SEDS UTILIZATION 

The weekly run of the CMI program provides CMI values for the 38 
stations shown in figure 3-1. These values are analyzed by the 
National Weather Service. The land analysis is then prepared 
for digitization, digitized, and input into SEDS as a magnetic 
tape, as describea in section 3.4. 
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TABLE .3-1. - SURFACE METEOROLOGICAL STATIONS 
USED IN DEVELOPMENT OF CMI 


STATION NAME 

STATION IRJMBER 

LATITUDE 

LONGITUDE 

ELEVATION 

ACAPULCO 

76805 

16 5 ON 

99 56W 

3 M 

CAMPECHE 

76695 

19 5IN 

90 32W 

25 K 

CHIHUAHUA 

76225 

28 38N 

106 04 w 

1423 M 

CIUDAD LAS CASAS 

7681;5 

16 44 n 

92 38W 

2128 'A 

CIUDAD LERDO 

76381 

25 SON 

103 32 W 

Il 40 M 

CIUDAD OBREGON 

767L6 

18 32 N 

92 39 W 

2 M 

COAIZACOALCOS 

76 Tkl 

18 09 N 

94 25 W 

l 4 M 

ENSENADA 

76050 

31 53 N 

116 38¥ 

15 M 

GUADALAJARA 

76612 

• 20 4 in 

103 20 W 

1589 M 

GUAITAJUATO 

76577 

21 OIN 

101 15W 

2073 M' 

GUAYMAS 

76255 

27 55 N 

110 53 W 

4 M 

JALAPA 

76687 

19 32 N 

96 55 W 

1399 M 

LA PAZ 

76U05 

24 ION 

110 12W 

18 M 

LEON 

76575 

21 OTN 

101 4 lW 

1809 M 

MANZANILLO 

7665L 

19 o 4 n 

104 20W 

8 M 

MAZATLAIJ 

76U58 

• 23 IIN 

106 25W 

78 M 

i-ERIDA 

76 Ghk 

20 58 N 

89 38W 

22 A 

MONCLOVA 

76342 

26 55N 

101 25W 

587 M 

MONTERREY 

76393 

25 40 N 

100 l 8 w 

534 M 

MORELIA 

76665 

19 42 N 

101 OTW 

1923 M 

OAXACA 

76775 

17 04 n 

96 43W 

1563 M 

PIEDRAS NEGRAS 

76243 

28 42 N 

100 32 W 

220 M 

PROGRESO 

76593 

21 ITN 

89 40 W 

iL M 

PUEBLA 

76685 

'19 02 N 

98 IIW 

2150 

SALINA CRUZ 

76833 

16 12N 

95 12 W 

56 M 

SALTILLO 

76390 

25 26n 

101 OOW 

1609 M 

TACUBAYA 

76679 

19 24 n 

99 11 ¥ 

2309 X 

TAvIPICO 

76548 

22 12 N 

97 51 W 

20 A 

TAPACHULA 

76903 

14 54 n 

92 l 6 w 

168 M 

TORHEON 

76382 

25 34 n 

103 25 W 

1150 

VERACRUZ 

76692 

19 12 N 

96 08w 

16 :• 

ZACATECAS 

76525 

22 4 th 

102 34 w 

2612 M 

BROWNSVILLE 

72250 

25 54 n 

97 26¥ 

10 

EL PASO 

72270 

31 48 n 

106 24 ¥ 

1200 M 

SAN DIEGO 

72290 

32 44 N 

IIT low 

11 M 

YUMA 

72280 

32 40 N 

ll 4 36 w 

62 M 


5-3 














The digital CMI values are used on a point-by-point basis in the 
SEDS data base to calculate a growth potential image and to 
prepare the combined product as described in section 2. 4. 1.1. 


3.2 DATA HANDLING 

The basic source o£ input data is the daily SEDS Surface 
Observations shown in table 2- IV. The data is summarized for 
Sunday through Saturday each week in the form of average tem- 
perature and total precipitation. The data is keypunched and 
submitted to the JSC Univac 1110 for processing. 

Detailed documentation of the data handling procedures may be 
obtained from C. M. Barnes, Ph. D. , D.V.M. 

3.3 VERIFICATION 

Verification of the CMI product depends upon indirect methods; 
of primary importance is feedback from field observations. 

Table 3- II giving a generalized’ interpretation of the CMI in 
terms of field conditions can be invaluable in both original 
analysis and verification. 

The CMI values change relatively slowly. Therefore, continuity 
in the CMI patterns should be expected. Discontinuity from the ' 
previous weeks analyses should be expected immediately following 
major occurrences of rainfall. 

The U.S. analysis of CMI is useful in the Mexican-U.S. border 
regions. The U.S. operational analysis is calculated from 
meteorological data that has been averaged over climatological 
divisions, representing area averaged CMI values as opposed to 
the SEDS point values. However, in the regions of overlap 
between the two analyses good agreement in analysis should be 
expected. 
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TABLE 3-II. - INTERPRETATION OF CMI VALUES 
IN TERMS OF FIELD CONDITIONS 


Index 

Interpretation 

Above 3 . 0 

Excessively wet, some flooding 

2.0 to 3.0 

Too wet, some standing water 

1.0 to 2.0 

Wet in spots 

0 to 1.0 

Normal 

0 to -1.0 

Topsoil moisture short 

-1.0 to -2.0 

Abnormally dry 

-2.0 to -3.0 

Too dry 

Below -3.0 

Exceptionally dry 


3.4 PREPARATION OF CMI TAPES 

There were several candidate methods for preparing digital 
tapes containing CMI data. All were designed to take raw CMI 
data for certain stations in Mexico calculated as described in 
section 3.2 to produce a tape containing CMI values for every 
pixel. 

The first alternative was the most general. It involved the 
interpretation of these data by a meteorological technician, 
who prepared an isoline map, which assigned CMI values over the 
entire SEDS grid. This image was then colored, read into a 
digital computer, the colored zones classified by standard 
techniques, and each pixel assigned a CMI value according to 
its class. This procedure is documented extensively in 
L. Giddings and E. Thompson, "An Interactive Method for 
Digitizing Zone Maps," CLEC-6498; JSC-09809^ . 
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This method was chosen for development since it was the most 
general. It allowed the preparation of digital tapes of any 
data that could be represented as zonal maps. To eliminate 
some dependence on a trained operator this method was largely 
automated, using standard LARSYS classification program. This 
required some modifications which are described in L. Giddings 
and S. Boston, "Digitizing Zone Maps using a Modified LARSYS 
Program," ChEC-7498; JSC-10757). "This modification converted 
long operation into a relatively short one. 

Other methods could also be used. The most feasible one would 
involve direct interpolations of the original data, or direct 
interpolation of data extracted from an interpreted image. For 
example, raw CMI values for selected stations (grid points) 
could be interpolated in two dimensions to fill the SEDS grid; 
alternatively, the raw values could be interpreted by a meteor- 
ologist, and data points extracted from the interpretation could 
be interpolated to fill the SEDS grid. Development of these 
methods has not yet been completed. 
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PROCESSING OF EMISSIVITY ZONE MAPS 


4.1 INTRODUCTION TO THE PROBLEM ‘ 

The SEDS screwworm algorithms are dependent upon the accuracy o£ 
the NOAA satellite measured ground surface temperatures. For 
accurate measurement, radiometric temperatures taken by the 
satellite must he corrected for certain errors, including the 
emissivity characteristics of the measured surface Cfig^^s 4-1). 
This section describes the emissivity correction utilized in 
SEDS and elaborates on the procedures for calculating the 
correction. 

The thermal infrared emissivity, e, of a homogeneous surface is 
defined (ref, 3) as the ratio of the emittance of the surface at 
a specified wavelength and emitting temperature to the emittance 
of an ideal blackbody at the same wavelength and temperature. 
Values for emissivity may range from zero to unity. 

To correct the problem, it is necessary to measure the emissiv- 
ity of each surface and use the measured value to correct any 
radiometric measurements of that surface. A complete analysis 
of the emissivity problem, and procedures for analysis are found 
in Phinney and Arp (ref. 4).. 

4.2 DEVELOPMENT OF THE EMISSIVITY ZONE MAP 

A map of emissivity values for every pixel of SEDS is necessary 
to provide the correct emissivity for every portion of Mexico. 
Since it is not possible to measure the emissivity of each por- 
tion of 'Mexico, it is necessary to have an efficient mapping 
technique . 

Analysis of the distribution of emissivity values and vegetation 
previously obtained in Veracruz reveals the following correlation. 
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TRUE 



ESTIMATED EMISSIVITY 

Figure 4-1. — Temperature error (°G) associated with an 
incorrect assumption, of emissivity at 300'^K. 



Variation in emissivity is primarily in response to vegetation 
and soils. In a study o£ emissivity by Phinney and Arp 
(ref. 5) , it was shown that visible and infrared data can be 
composited to produce a zonation of ground data. These zones 
proved to correspond to the distribution of vegetation in the 
Veracruz test site. By knowing the distribution of various 
plant communities (figure 4-2) and their calculated community 
emissivity values, it is possible to synthesize a map of the 
emissivity zones of Mexico, 

4.3 METHODS OF PRODUCTION OF THE SEDS EMISSIVITY MAP 

4.3.1 SOURCES OF EMISSIVITY MEASUREMENTS 

In accordance with procedures set forth in Phinney and Arp 
(ref. 4), emissivity data were field gathered from each major 
vegetation community and soil type in Mexico. Determinations of 
the important communities of plants in Mexico were provided by 
analysis of the works of G, Flores-Mata et al . (ref. 6); J. S. 
Beard (ref. 7 ^ 8); F. Miranda and E. Hernandez (ref. 9) 1963; 

A. S. Starker (ref. 10); and S. R. Eyr.e (ref. 11). 

4.3.2 SOURCES OF RADIOMETRIC DATA 

Each day NOAA sends two raw data tapes in analog form to NASA 
JSC. One tape contains the day visible and day infrared (IR) 
data; the other tape contains the night IR data. These tapes 
are registered and converted by Aeronutronic Ford from analog 
to SEDS digital form to produce three registered data tapes, day 
IR, day visible, and night IR. The three daily tapes are trans- 
ferred to the Health Applications Office in Universal format 
suitable for direct read on the Aerojet Data Analysis Station 
(DAS) . Day IR data are received at full resolution, while day 
visible and night IR data are compressed. At least two weeks of 
data are amassed to assure enough data to create a cloud-free 
image in all three channels. 
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Tapes are reformatted on the Aerojet DAS from Universal to 
LARSYS format using the standard system program U2L3AS and com- 
bined to make a multifile tape. 

4.3.3 SELECTION OF IMAGES FOR USE IN COMPOSITING 

The best images for use in compositing are selected by either 
screening each, image on the Aerojet DAS and photographing with a 
Polaroid camera, or by viewing lantern slides of the SEDS data 
prepared by Aeronutronic Ford. In either case, images are 
selected so that no two images have clouds in the same regions. 
In non-cloudy periods, as few as two images may be used to 
create a cloud-free composite. In cloudy periods, five raw 
images may be needed to composite each data channel. 

4.3.4 CREATION OF COMPOSITED IMAGES TAPE 

After the desired images are selected, they are composited onto 
a disk of the Aerojet DAS by use of the E-DISK routine. For 
program documentation, contact Dr. C. M. Barnes, Health Applica- 
tions Office. The visible data is composited using a minimum 
brightness test, whereby each pixel of the second image is com- 
posited with its complementary pixel on the first image. The 
pixel that has the minimum brightness is recorded, thereby 
deleting the brightest cloud-contaminated pixels. Infrared data 
is similarly composited on separate files of the disk. In this 
case, each pixel of the second image is compared with the com- 
plementary pixel of the first image using a maximum temperature 
criterion. In this case, the warmer of the two pixels is 
recorded. This process removes any pixels contaminated by cold 
clouds. Completion of the compositing yields a 3-channel disk 
which is immediately read onto a 9-track tape, yielding a 
3-channel, 9-track tape in LARSYS .format for analysis. 


4-5 



4.3.5 CLUSTER ANALYSIS OP 'COMPOSITED DATA 

The composited data tape is treated as a 3-channel multispectral 
scanner (MSS) tape. The data tape is analyzed using the ISOCLS 
processor on the UNIVAC 1108 to define those areas of similar 
temperature and reflectivity. The clustering procedure used in 
ISOCLS (ref. 12) is an iterative procedure, which assigns each 
data sample to a specific cluster by determining the nearest 
cluster center (in absolute distance) and assigning the sample 
to it . 

To begin clustering, all program defaults are used. For the 
radiometric data under study, about 20 classes will be created 
by the ISOCLS processor. Of the 20 classes, several will be 
water classes and several will be cloud classes. To determine 
the identity of the remaining classes, the interaction of a 
trained botanist is necessary. Frequently, two or more computer 
defined classes will correspond to various phases of a single 
vegetational zone. On occasion, the computer will not distin- 
quish separate vegetative zones. This is frequently due to 
seasonality effects, but careful study is necessary to avoid 
serious misclassif ication of the data. 

This process results in the creation of several clusters or 
regions which have similar temperature and reflectivity values. 
These regions correspond to vegetation zones and are ultimately 
presented in a map-lihe format to which emissivity values are 
assigned. The clustered results return as a statistics deck 
and 110 lines of graymap output. Initial verification of the 
clustered data based upon ground truth evidence may proceed at 
this point using the graymapped output. However, to receive a 
complete output of the classified image, the composited data 
must be processed again. 
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4.3.6 PRODUCTION OF THE FIRST HARD COPY OF THE EMISSIVITY ZONES 


The ISOCLS processor only produces a minimal 110 lines o£ output 
for use in interpretation of the clustered data and only analy- 
zes every second or fourth element for purposes of efficiency. 

To receive all the clustered data, the composited data tape 
(9-track) plus the previously computed statistics deck are run 
on the UNIVAC 1108 through the CLASSIFY processor (ref. 12), 
which classifies the entire image data based upon. the statistics 
computed form the ISOCLS processor. This procedure produces a 
7-track tape "MAPTAP" . containing every pixel classified accord- 
ing to the statistics generated by ISOCLS. 

4.3.7 PRODUCT DISPLAY 

The 7-track MAPTAP is run through the 1108 again to produce the 
first hard copy of the entire image by using a program called 
DISPLAY. Results of the DISPLAY routine include a 9-track DAS 
tape for use on the Aerojet DAS, in addition to the printouts of 
the entire study area (ref. 12). 


4.3.8 ASSIGNMENT OF EMISSIVITY VALUES TO ZONES 

The DAS tape (9-track) product from the DISPLAY routine is dis- 
played on the Aerojet DAS using the Color Study program. Each 
zone defined in the ISOCLS routine is displayed and analyzed 
against known ground truth information. In some cases, a 
cluster represents a discrete zone on the ground; in other 
cases, it comprises only part of a larger zone. Identification 
of zones made from field work and references, and computation 
of the community emis.sivity for each zone are described in 
Phinney and Arp (ref. 4) . Assignment of emissivity values is 
made to each zone by use of the E-MISS program using the formula 


f(e) = Integer 
e ave = .975 


( e ave 
e true 



(853.3333) + 80.5 
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Numbers for each community and the average emissivity value are 
encoded. The 9-tracTc LARSYS tape is then converted to Universal 
format and transferred to Aeronutronic Ford for use in processing 
of the daily radiometric data. For program documentation, con- 
tact Dr. C. M. Barnes, Health Applications Group. 

4.3.9 CURRENT STATUS 

Two tapes have been created to date. The spring tape uses 
selected data from 15 to 29 March and the fall tape uses selected 
data from 28 September to 9 October. In theory, no new emissiv- 
ity tapes need to be made since the vegetation of Mexico will not 
change measureably in the next year or two. 

In practice, however, it is not certain whether or not new tapes 
will need to be made. Very wet seasons or extended droughts will- 
affect the vegetation strongly. There is no guarantee that the 
current tapes are typical for Mexico, or that this year was 
normal in all respects. In all probability, at least two more 
tapes should be generated for 1976 to determine ^^^hat changes, if 
any, may result from year to year. 
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